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ABSTRACT 

The unanticipated discovery of the first close-in planet around 51 Peg has 
rekindled the notion that shortly after their formation outside the snow line, some 
planets may have migrated to the proximity of their host stars because of their 
tidal interaction with their nascent disks. After a decade of discoveries, nearly 20 
% of the 200 known planets have similar short periods. If these planets indeed 
migrated to their present-day location, their survival would require a halting 
mechanism in the proximity of their host stars. Here we consider the possibility 
that a magnetic coupling between young stars and planets could quench the 
planet's orbital evolution. Most T Tauri stars have magnetic fields of several 
thousand gausses on their surface which can clear out a cavity in the innermost 
regions of their circumstellar disks and impose magnetic induction on the nearby 
young planets. After a brief discussion of the complexity of the full problem, we 
focus our discussion on evaluating the permeation and ohmic dissipation of the 
time dependent component of the stellar magnetic field in the planet's interior. 
Adopting a model first introduced by Campbell for interacting binary stars, we 
determine the modulation of the planetary response to the tilted magnetic field 
of a non-synchronously spinning star. We first compute the conductivity in the 
young planets, which indicates that the stellar field can penetrate well into the 
planet's envelope in a synodic period. For various orbital configurations, we 
show that the energy dissipation rate inside the planet is sufficient to induce 
short-period planets to inflate. This process results in mass loss via Roche lobe 
overflow and in the halting of the planet's orbital migration. 
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1. Introduction 

Perhaps the most surprising finding in the search for extrasolar planets is the discovery of 
short-period (P < 1 week) Jupiter-mass (Mj) companions around Solar-type main sequence 
stars (Mayor & Queloz 1995, Marcy et al. 2000). Among the inventory of >200 presently- 
known extrasolar planets, 20% have P = 1 — 7 days. Nearly 20 short-period planets have 
measured radii (R p ) that are comparable to or larger than that of Jupiter (Rj). While these 
information may be biazed because of observational selection effects, these planets are most 
probably gas giants. 

According to the conventional sequential-accretion scenario (Pollack et al. 1996), the 
most likely birth place for gas giant planets is just outside the snow line where volatile heavy 
elements can condense and coagulate into large planet building blocks (Ida & Lin 2004). In 
protostellar disks with comparable surface density (£), metallicity ([Fe/H]), and temperature 
(T) distributions as those of minimum mass nebula model (Hayashi et al. 1985), protoplanets 
with M p ~ Mj induce the formation of a gap near their orbit as a consequence of their tidal 
torque on the nascent disks (Goldreich & Tremaine 1978, 1980, Lin & Papaloizou 1980, 
1986a, 1993). In relatively massive and fast evolving disks, the outward transfer of angular 
momentum due to the disks' intrinsic turbulence can lead to an inward mass flux (M d ) and 
the migration of the gas giant planets (Lin & Papaloizou 1986b). This process is commonly 
referred to as type II migration (Ward 1997). 

This migration scenario was resurrected to account for the origin of the first known short- 
period extra solar planet (Lin et al. 1996). Although type-II migration provides a natural 
avenue for relocating some gas giants, a mechanism is needed to retain these planets close to 
their host stars. Moreover, many stars are born with rapid rotation (Stassun 2001). When 
young planets venture close to their host stars, angular momentum would be transferred from 
the stellar spin to the planet's orbit if the stellar spin frequency is still larger than the 
planet's orbital frequency (f2fc). The rate of the star-to-planet angular momentum transfer 
intensifies rapidly and may exceed that from the planet to the disk. 

Two basic physical effects were suggested as potential migration barriers. The first one is 
tidal interaction between the host star and the planet. The gravitational perturbation of the 
star and close-in planet leads to responses in both the star and planet. For tidal frequencies 
smaller than twice the spin frequency, inertial waves are excited in the convective envelope 
of the host star and are dissipated there by turbulent viscosity (Ogilvie & Lin 2007). But, 
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the tide excited by a close-in gas giant planet in a star, with a structure similar to that of 
the present Sun, marginally fails to achieve nonlinearity so that their survival is ensured. 
Nevertheless, during the formation epoch of solar-type stars, conditions at the center of the 
star evolve, so that nonlinearity may set in at a critical age, resulting in a relatively intense 
star-planet tidal interaction. 

The second effect suggested is based on the magnetic interaction between the host star 
and the planet. Young T Tauri stars also have radii (R*) 2-3 times that of the present-day 
Sun (Rq) and several thousand gausses fields (B*) on their surface (Johns-Krull 2007). The 
stellar magnetosphere threads across the inner regions of the disk and clears a cavity out 
to a critical radius (R c ) which is determined by both the magnitude B* and Md (Konigl 
1991). The subsequent complex interplay between accretion and outflow leads to angular 
momentum exchange which induces u;* to evolve toward Q k at R c (Shu 1994). When the 
planet's orbital semi major axis (a) reduces well inside R c , its Lindblad resonances relocate 
inside the star's magnetospheric cavity. In principle, the planet's migration would stall due 
to its diminishing tidal torque on the disk. 

However, if the star's magnetospheric interaction with the disk can lead to = Q k (R c ), 
the planet inside the magnetospheric cavity would have Q k > In this limit, the star-planet 
tidal interaction would induce a transfer of angular momentum from the planet to the star. In 
addition, the differential motion between the planet and the stellar spinning magnetosphere 
induces an electromagnetic field with a potential to generate a large current analogeous to 
the interaction between the Jovian magnetosphere with its satellite Io (Goldreich & Lynden- 
Bell 1969). The associated Lorentz force drives an orbital evolution toward a synchroneous 
state, in which case, angular momentum would be transferred from the planets with Q k > ^* 
to their host stars, and the planets would continue their orbital decay. 

In order to determine the necessary condition for the retention of close-in young plan- 
ets, we examine, in this paper, their interaction with the magnetosphere of their host T 
Tauri stars. In §2, we briefly recapitulate the essential concepts and validity of previous 
investigations on some related topics and give an overview of the key phenomena that will 
be discussed in this and later papers. In §3, we adopt an existing model in order to ex- 
amine the interaction between a planet and the magnetosphere of its host star. In §4, we 
compute precisely the planet's magnetic diffusivity for a specific set of parameters, as well 
as the corresponding ohmic dissipation rate within that planet. In §5, we suggest that the 
ohmic dissipation can generate sufficient heat to inflat the planet. In §6, we construct an 
idealized self-consistent model in which the polytropic and isothermal equations of state are 
utilized. These equations represent the expected outcome of radiation transfer within the 
fully convective interior and the isothermal surface of a close-in planet which is exposed to 
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the intense radiation from its host star. This internal model allows us to compute the mag- 
netic diffusivity. With these tools, we discuss in §7 the structural adjustment of the planet 
in response to this heating source and we compute the ohmic dissipation and mass loss rate 
for different set of parameters. Finally in §8, we summarize our results and discuss their 
implications. 

2. Planetary and astrophysical analogue 

Two previous analyses are directly relevant to the present study: 1) the interaction of 
Io with the magnetosphere of Jupiter and 2) the spin-orbit synchronization in binary stars 
containing a magnetized white dwarf and its main sequence or white dwarf or planetary 
companion. 

2.1. Unipolar Induction in Io 

Io orbits around Jupiter inside its magnetosphere once every 1.7 days which is consid- 
erably longer than Jupiter's 10 hours spin period. This relative motion imposes a periodic 
variation in Jupiter's decametric emission (Duncan 1966). A class of models that accounts 
for the origin of this emission was developped based on the assumption that Io has a suf- 
ficiently high conductivity. In Io's rest frame, the electric field vanishes and the steady 
component of Jupiter's magnetic field permeates in Io's interior over time. When a steady 
state is established, the tube of constant magnetic flux is firmly frozen into Io (Piddington 
and Drake 1968) due to its high conductivity. The flux tube carried by Io moves through 
the surrounding field lines (which corotate with Jupiter) and slips through Jupiter's less 
conductive ionosopheric surface. Plasma in Jupiter's ionosphere flows around the tube and 
introduces a potential difference across it (Goldreich & Lynden-Bell 1969). The associated 
electric field drives a current which travels down one half of the flux tube from Io and is sent 
back to Io along the other half. Within the flux tube connecting Io and Jupiter's ionosphere 
and those across it on Io, the electric field vanishes as a consequence of high conductivity. 
Thus, this DC circuit is closed by Io as a unipolar inductor. 

The magnitude of the electric current is primarily determined by the Pederson conduc- 
tivity at the foot of the flux tube i.e. on the ionosphere of Jupiter. Finite conductivity also 
determines the magnitude of the drag against the slipage of the flux tube through Jupiter. 
This drag results in energy dissipation on the surface of Jupiter and in a torque on the 
orbital motion of Io, driving the system towards a state of synchronization. This configu- 
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ration is justified by the assumption that a constant flux tube is firmly anchored on and 
dragged along by Io which requires the conductivity in Io to be much higher than that in 
Jupiter's ionosphere. A 10° inclination between Jupiter's magnetic dipole and rotation axis 
does introduce a periodic variation (over a synodic period) in the field felt by Io. The per- 
meation and dissipation of this time-dependent AC field may be negligible in the limit of 
high conductivity in Io. 

The validity of the key assumption for the unipolar induction model (i.e. conductivity 
on Io is larger than that in Jupiter's ionosphere) has also been challenged by Dermott (1970). 
A modest resistence in Io would distort the field which may lead to field slippage through 
Io. In this case, the passage of Io through the magnetosphere of Jupiter would lead to the 
generation of Alfven waves along the flux tube (Drell et al. 1965, Neubauer 1980). But, 
due to the field displacement, the waves, partially reflected at the foot of the flux tube on 
Jupiter's surface, may not be able to return to Io, in which case the DC circuit would be 
broken and the motion of Io would be decoupled from the that of the flux tube. Nevertheless, 
the Alfven waves are dissipated inside both Io and Jupiter, leading to a torque which must 
depend on their penetration depth. 

An alternative class of scenarios has been proposed based on the assumption that the 
magnetosphere is everywhere anchored on Jupiter and the flux tube moves freely through Io 
(Gurnett 1972). This model requires the conductivity in Jupiter's ionosphere to be larger 
than that in Io. It assumes that the presence of Io creates a plasma sheath with an electric 
field to cancel the induced EMF associated with the motion of Io relative to Jupiter's magne- 
tosphere (Shawhan 1976). The simplifying approximations in the development of this theory 
have been challenged by Piddington (1977) who questioned both the validity of the sheath- 
creation mechanism and the self-consistency of the internal and external field configurations, 
subjected to the electric currents in and around Io. 

On the observational side, UV emissions from Io's footprint on Jupiter has been ob- 
served. But, it extends well beyond the intersection between Io's flux tube and Jupiter's 
ionosphere and the emssion downstream is protracted (Clarke et al. 1996). These obser- 
vations do not agree with the simple interpretation of either the unipolar induction or the 
plasma sheath scenarios. 

2.2. Magnetic Coupling in Interacting Binary Stars 

There are many close binary-star systems with a white dwarf as their primary compo- 
nent. These systems also contain main sequence stars and other white dwarfs as secondary 
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components in compact and circular orbits around each other. In some cases, mass is trans- 
ferred from the secondary to the primary. In other cases, gravitational radiation may play 
an important role in determining the evolution of these systems. 

A sub class of such interacting binary stars, AM Her systems, is composed of a magne- 
tized white-dwarf primary and a lower-mass main sequence star as its secondary in a fully 
synchroneous orbit despite the ongoing mass transfer between them (Warner 1995). This 
orbital configuration is very similar to that of the Jupiter-Io system despite the enormous 
difference between the mass ratio in the two cases. The motivation for studying the impact 
of magnetic coupling between these stellar components is to assess whether this synchronous 
state can be achieved through the ohmic dissipation of the white dwarf's field in the main 
sequence star's surface (Joss et al. 1979). Toward this goal, Campbell (1983, hereafter C83) 
adopted a novel approach by considering the penetration and dissipation of a periodically 
variable field, associated with an asynchronously spinning primary. 

Campbell's approach is fundamentally different from that of the unipolar induction 
model. In this analysis, Campbell focussed on the flow in the envelope of the secondary 
and neglected the possibility of current flowing through the flux tube between the secondary 
star/satellite and the surface of the primary star/planet. This vacuum-surrounding approxi- 
mation is justifiable since conductivity on the primary is likely to be much larger than that on 
the secondary and the stationary component of the field is frozen in the white-dwarf primary 
but not in the main-sequence secondary. Campbell analyzed the time dependent response 
of the secondary, including the modification of the field by the induced (AC) current in it 
(Campbell 2005), to the periodic modulation of the field. In contrast, the unipolar induc- 
tion model depends on the explicit assumption that the field is anchored on the secondary 
and its distortion near the secondary must be small so that a complete current loop can 
be established between the primary and the secondary. Campbell determined the periodic 
diffusion of the field and the ohmic dissipation of the induced AC current in the companion 
whereas that of the induced DC current is assumed to occur in the primary in the unipolar 
induction model. 

In recent applications of the unipolar induction model in the context of interaction 
between white dwarf binary stars, the current induced by the unperturbed field has been 
computed but the induced field generated by the current was neglected (Wu et al. 2002, 
Dall'Osso et al. 2006). A tidal torque is computed at the footprint of the flux tube, which is 
attached to the secondary white dwarf, on the surface of the magnetized primary white dwarf. 
A totally self-consistent solution of this diffult and complex problem remains outstanding. In 
addition to the uncertain anchorage location of the field, it is not clear whether the resulting 
misalignment of the total (original plus induced) field and current may be sufficiently large 
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to break the circuit, in which case, Campbell's model may be more appropriate. 

2.3. Mathematical approximations made by the two previous models 

In this subsection, we summarize the physical description of the two models just pre- 
sented (the unipolar inductor versus the periodic diffusion) by explicating the mathematical 
approximations made in each one of them. The complete MHD induction equation can be 
expressed as 



where the magnetic diffusivity 77 = l/fj,o<7 and the electrical conductivity a functions of 
position and /xo is the permeability. If rj is constant, the second term on the right hand 
side becomes —V A (77V A B) = 77 V 2 B which reduces to a common expression of diffusion. 
The two models (unipolar induction versus periodic diffusion) consider two complementary 
approximations of equation [TJ In the problem where Io is treated as an unipolar inductor, 
its conductivity is explicitly assumed to be large so that the second term on the right hand 
(the diffusion term) is negligible compared to the first (i.e. the induction term). In this 
configuration, one can show that the field lines of the steady component of the magnetic 
field are moving with Io and appear to be "frozen" on Io (see appendix A). Alternatively, in 
the model considered by Campbell, it is the first term on the right hand side that is being 
neglected. Moreover, only the diffusion of the time dependent component of the field is being 
considered. This approximation is valid if the two interacting bodies are almost in corotation 
(i.e. the relative speed v that appears in equation (TTJ is small), or if the conductivity in the 
secondary is small. 



The process under investigation in this paper is analogeous to both the Jupiter-Io and 
the interacting binary star problems. In fact, the unipolar induction model has already been 
applied to study the orbital evolution of terrestrial planets or the cores of gas giants around 
white dwarfs (Li et al. 1998). There are even follow-up determination of the radio flux 
densities from potential white-dwarf /planet systems (Willes & Wu 2005). In this analysis, 
although the dissipation of the induced current due to the finite conductivities in the white 
dwarf was considered, the feedback modification of the field and the dissipation within the 
planet have been neglected (Li et al. 1998). As discussed above, it is not clear whether a 
DC circuit can be closed to promote the unipolar induction mechanism. 





2.4. Overview of the phenomena that will be discussed 
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In light of these uncertainties, we consider both classes of models for the interaction of 
close-in planets with their magnetized host stars. In this paper, we focus our discussion on 
the mechanism described by Campbell, and apply it to a hot-jupiter revolving around its 
star. We will return to the unipolar induction problem in a later paper. 

When young planets first arrive at the vicinity of their host stars, they are unlikely 
to be in a totally synchronized state. The stellar magnetic field felt by the planet may 
be dominated by the periodic modulation associated with the synodic (between the stellar 
spin and the planet's orbit) motion. In addition, the temperature in the planet's surface is 
expected to be ~ 10 3 K and the conductivity there may be moderate. In response to the 
modulation of the field, the interior of the planet continually adjusts to the magnetization 
effects so that the flux tube cannot be effectively frozen in the planet. All of these boundary 
condition suggest that at least over some regions of the planet (especially on the night side 
where the photo ionization due to the stellar flux is negligible), the modulation of the field 
may lead to an induced current inside the planet which does not contribute to the close 
circuit of a unipolar inductor. 

Following the geometry introduced by Campbell (C83), we consider a close-in gas giant 
planet, with a finite conductivity, interacting with a time-dependent magnetic field generated 
by the star. An induced current is generated inside the planet, which is associated with an 
ohmic dissipation rate. Our main contributions to the model used by Campbell are: 1) the 
relevant diffusivity inside the gas giant planets, 2) the effects of the ohmic dissipation on the 
planet's internal structure, and 3) the resulting orbital evolution of the planet. (Items 2 and 3 
have negligible consequence in the interacting binary star problem considered by Campbell). 
Since we are only considering the dissipation in the planet's interior, the associated torque 
applied on its orbit should be regarded as a lower limit. 

In our scenario, we postulate that at sufficiently close proximity to the host star, the 
stellar magnetic field is sufficiently intense that the ohmic dissipation of the periodically 
diffused field inside the planet is adequate to heat and to inflat the planet until it overflows 
its Roche lobe. The hemisphere of the planet facing its host star is also exposed to the 
intense flux of UV radiation during the stellar infancy. It is possible for the planet to 
develop a substantial ionosphere regardless the state of synchronization between the planet's 
orbit and spin (the time scale for establishing local ionization equilibrium is much faster 
than the planet's spin and orbital periods). 

We will separately study these two phenomena (angular momentum transfert due to 
mass loss and presence of a ionosphere) in the follow-up papers of this serie. We will show 
that the angular momentum transfer associated with the mass transfer can halt the orbital 
evolution of the planet. We will also present an analysis on the conductivities in the planet's 
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day-side ionosphere and on the host stars surface. This will lead to an analysis on the 
condition for the unipolar induction to effectively operate and apply a significant slow-down 
torque on the planet's orbit. 

3. Magnetic induction 

In this section, we are going to derive the governing equations that we use to compute the 
ohmic dissipation rate. Various equations are presented here for the purpose of introducing 
the algorithm of the numerical models to be presented in subsequent chapters. Although 
we follow closely the approach made in C83, for brevity, we do not repeatedly cite this 
reference. But, wherever similarities occur, referal of Campbell's earlier work is implicitly 
implied. Also, throughout the paper, we use SI units. 




Fig. 1. — The geometry of the system. The star is on the left, at the center of the set of 
axies (xo,yo,zo), and the planet is on the right, at the center of the set of axies (x,y,z). 

We consider a protoplanetary system with a gas giant planet revolving around a T Tauri 
star with an angular frequency Q p . Well beyond the planet's semi major axis, there is also 
a protoplanetary disk. The host star has a dipolar moment m tilted with an angle a with 
respect to its spinning axis (see figure 1). The angular frequency of the stellar spin is a;*. The 
orbital axis of the disk and the planet are parallel to the star's spinning axis. The following 
analysis is applied to a frame of reference centered on the star and rotating with the planet. 

In this frame, the planet is a fixed object (the planet's spin is neglected) in a periodic 
magnetic field with a frequency uj = uu* — Q p . From Ohm's law J = crE and Maxwell's 
equations = —V A E and V A B = /i J, the equation on the magnetic field becomes: 

_ = -VA(7?VAB). (2) 
It follows that in the mechanism considered by Campbell (as well as in this paper), it is 
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the time dependent stellar magnetic field, diffusing inside the secondary (for Campbell) or 
the hot Jupiter (in our paper), as well as the planet's induced magnetic field, that generate 
the current inside the planet, following the equation V A B = fi 3. The relative speed 
between the planet and the stellar magnetic field thus intervenes not through E = — v A B 
but through the time dependence in the stellar magnetic field that diffuses in the planet. 

Following C83, we only consider the poloidal component (ft of the magnetic field: 

B = VA(VA(0e r )) (3) 

where (ft is a function of r, 6, if. and can be expanded in terms of the sperical harmonics 
Y™(9, (p) (equation H]). Moreover the variation in time of the magnetic field felt by the planet 
is periodic. In the limit where the field penetrates quickly in the planet compared to the 
time scale on which the field changes (so that the planet can respond "adiabatically" ) , we 
can account for the time dependence of (ft by mutiplying its spatial part by e iuJt : 



[r,t) = no 



CrG l {r)Yr{e^) 



Lm 



e luJt (1^0 and - I < m < /) (4) 



where C z m are constant coefficients and Gi(r) is a function of r to be determined. We then 
replace B in the left hand of (2) by its expression in (3). After integration, we obtain 

V AB = — — V A {(ft e r ). (5) 
V 

We then replace B in the left hand side of this equation using (jHJ), and develop both sides 
of the equation. After identification, we obtain: 

-rr(r) - v - ' H GAr) = inside the planet. 6 

dr 1 \ r 2 n J 

This equation holds inside and outside the planet (same as eqs 16 and 18 in C83). However, 
outside the planet, the conductivity is assumed to be very low and, therefore, the magnetic 
diffusivity is extremely high compared to the diffusivity inside the planet. In the limit where 
the magnetic diffusivity outside tends to infinity (equivalent to a vacuum surrounding), the 
equation (jSJ) becomes: 

^nr(r) - ( ^ GAr) = outside the planet. (7) 

dr 2 \ r 2 I 



We consider the radial part of the poloidal scalar outside the planet. Following C83, 
we introduce (ft s tar the radial part of the poloidal scalar outside the planet due to the star's 
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magnetic field, and 4> p ianet the radial part of the poloidal scalar outside the planet due to the 
planet (cf C83 eqs 21-22): 



'star 



u msma 2 . . . . * u msma , r n . 
— r [Zcosif smut + sm^ coscut) P, H — r PI ismut 



8vr d A 

— ^-cos2cp sincut + -sin2<yj coscut^ P 2 2 



(8) 



Jplanet 



+ 



(Pi sincut + p 2 coscut) + iaqP% 



cosy 9 / • \ smif . . 
[a\ smut + a 2 coscut) H (a 3 smut + «4 coscut) 

cos2cp 



/ • \ sin2cp . 
,, (71 sincut + 72 coscut) H — ^sincut + 74 coscut) (9) 



where = -sin0, P 2 ° = ±(3cos 2 fl 



1), and P| 



3sin 2 # are the associated Legendre 



polynomials (our convention for P/ has an opposite sign as that adopted by Campbell). In 
addition, 4> p i a net has the same time and angular dependence as 4> s tar because the field inside 
the planet is induced by the stellar 's magnetic field. 

The sum 4> ou tside = (ft star + (fipianet is the total poloidal scalar outside the planet, and 
4> outside (given by (jHJ) and Qj) is equal to (fiinside (given by (j4])) at the surface of the planet 
(r = PJ. 



3.1. Poloidal scalar inside the planet 

In order to determine the poloidal scalar inside the planet, we first numerically calculate 
the values of Gi(r) (the radial part of 0, cf. equation H]) and G[{r) inside the planet by solving 
equation (jfJJ). We then calculate the coefficients C™, which appears in the decomposition 
of </>. They are determined by the boundary conditions which connect the interior and 
exterior solutions. In the rest of this section (§3), we assume that the conductivity profile 
is known, and we describe the procedure used to compute the ohmic dissipation rate inside 
the planet. In the following sections, we apply the method described in §3 to compute the 
ohmic dissipation rate inside the planet. 

In the following sections, we compute the magnetic dissipation and numerically obtain 
the value of the ohmic dissipation rate. 
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3.1.1. Computation of G(r) 

If the diffusivity r)(r) is known, we can solve ([6]) numerically, for I = 1 and I = 2, 
with a two-point boundary solver using the Newton-Raphson-Kantorovich method, and the 
equations and boundary conditions are given below: 



Y x (t) = Ke(G{r)) 
Y 2 (r) = lm(G(r)) 



Equations : 

Y{(r) = Y 3 (r) 
F 2 '(r) = y 4 (r) 



F 3 '(r 



Yl(r) = ^(r) + 



Yo(r) 



Boundary conditions at r = R p & 
G{(i2p) + - (2/ + l)R l p = (10) 

G'(r ~ 0) - l -^Gi(r ~ 0) = 



Computation of the C(l,m) 

The complex coefficients CJ™ = /i™ + ii/J™ have real and imaginary parts [if 1 = Re(Cl n ) 
and v™ = ImiC™). We equate the real part of the decomposition of the poloidal scalar inside 
the planet given in (j3j) at r = R p (radius of the planet) with the expression of 4> ou tside = 
4>star + (ppianet given in flSJ) and <Q ai r = R p . 

Moreover, using the fact that (P^P^P^) and then (cosutcosip, coswisinyj, sinwicos^, 
smutsimp) (cosu;£,sinu;£), (cosu;£cos2<^, cosa;tsin2(^,sina;tcos2(^, smutsm2ip) are a set of bases, 
we get a set of linear equations which can be solved for fi^ 1 , v\, v x , an, a 2 ,a 3 ,a 4 ), 
v\, Pijfa), and /i^~ 2 , 7i, 72, 73, 74) (the linear systems verified by these unknowns 

are given in appendix B). 



3.2. Computation of the ohmic energy dissipation rate 

The potential generates an electric field E which induces a volumic current J inside the 
planet. The associated ohmic dissipation inside the planet is V vo iumic = -Re (J) -Re(E). Using 
E = -J and J = —V A B, we can write: 



V 



a r 



(TZe(3)) 2 dV 



a(r) 



Tie 



V AB 



.dV 



Moreover, using equation (jHJ), we can write 



P = — 



1 / dlm((j)) 



sin# \ dcp 



r, 9, ip) 



sm8 



( dlm(<f>) x 2 



dr d9 dip. (12) 



13 



We use eq. fl4]) to express the real and imaginary parts of $. After integrating over 9 and 
(p, we are left with: 

V = I < V volumic > r 2 dr (13) 



where the angle-integrated volumic power and 



^ r volumic ^ — ~ i COS tot 



(A 2 n + A\ 2 ) + 3 {A\ 7 + A 



' A 2 

.-"■15 



-sm 2 ujt 



(^13 + Au) + 3 (-A19 + Ajo) H A^g 



+sinwteoso;t 



(A12A14 + AuAxa) + 3 (A17A19 + A 18 A 2 



-A 15 A 



16 



where the expressions for Aj,- are given in Appendix C. 



4. Conductivity profile and ohmic dissipation rate 

The general setting of the problem and the basic equations have been laid down. We 
have seen that once a conductivity profile is chosen, one can solve (TTU|) and determine Gi(r) 
inside the planet (the radial part of the poloidal scalar inside the planet). Then, one can 
compute the C(l,m), and finally obtain the ohmic dissipation rate V inside the planet. 



4.1. Computation of V for one specific set of parameters 

We compute the conductivity inside the planet with two parallel approaches. In §§5-7, 
we develop an idealized self-consistent internal structure model to determine the response of 
the planet to the ohmic dissipation of the induced current in it. But, in this section, we first 
introduce a realistic, but non self consistent, model with the following set of parameters: 



Planet's mass and radius: 0.63M/ & R p = 1.4 Rj = 10 8 m, 
Semi-major axis: a = 0.04 AU = 6 x 10 9 m. 



These and other (such as mass and luminosity) stellar parameters are appropriate for the 
short-period planet around HD209458 (Bodenheimer et al. 2001). We compute the internal 
conductivity due to the ionization of the alkaline metals (see Appendix D for details). Al- 
though the planet is heated on the day side, thermal circulation can redistribute the heat 
and reduce the temperature gradient between the two side of the planet (Burkert et al. 2005, 
Dobbs-Dixon & Lin 2007). We adopt a spherically symmetric approximation for the surface 
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temperature of the planet to be 1,360K. Here, we neglect the modification in the internal 
structure due to the ohmic dissipation which will be considered with self-consistent models 
in the next sections. In Paper IV, we will also consider the conductivity on the planet's 
upper atmosphere due to photoionization which only occurs on the dayside of the planet. 

Using this conductivity profile, we can approximate the magnetic diffusivity i](r) = 
l//i <r(r) by: 

r]{r) ~ 10 3 exp 

where the effects of the photoionization have been neglected in this paper. 
To apply the procedure described in §3, we also need to specify: 

Relative angular velocity: uj = 10 5 s _1 , 

Star's magnetic dipole: m = 4 x 10 34 A m 2 , 

Value of the tilt of the magnetic dipole: sin (a) = I. 

We then obtain the following V (also see figure 2 for the graphs of Gi(r)) 

V(t) = 2.26 x 10 21 cos 2 cut + 2.1 x 10 21 sin 2 cut + 1.3 x 10 21 sinwtcosu;t (15) 
We then take the average in time over one synodic period and obtain V ~ 2.18 x 10 21 Js^ 1 

The conductivity profile we have obtained here is sensitive to the planetary structure 
model. At the epoch of planet formation, the gas accretion and planetesimal bombardment 
history are stochastic (Zhou & Lin 2007). The opacity in the accretion envelope of proto gas 
giant planets may also be subjected to variations due to dust coagulation (Iaroslavitz et al. 
2007). The thermal evolution of these planets can be highly diverse. There may, therefore, 
be a dispersion in the magnitude of rj. 

4.2. Comments on the skin depth and the dependence of the ohmic 
dissipation on the conductivity and on the sign of u 

Once the conductivity profile within the planet is determined, we are able to compute the 
energy dissipation rate inside the planet of the current induced by the star's magnetic field. 
In light of the possible uncertainties in the magnitude of rj, we compute the ohmic dissipation 
rate for different rj by artificially modifying the above determined rj with a multiplicative 
factor. The resulting magnitude of the time-averaged value of V is listed below (tabled]). 
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Im(Gi = 2){r), and their first derivatives, 

2" 



10 d exp 



25 



The shape of Gi for 



I = 1 and I = 2 are very close but the amplitudes for I = 2 are about 10 8 = R p higher than 
for I = 1. Indeed, the major difference between I = 1 and I = 2 is found in the equations 
describing the boundary conduitions (see equation [TUI where a factor 10 8 between I = 1 and 
I = 2 comes from the term R l ). In addition, we found \C(l,m)\ for I = 1 is about 10 10 times 
larger than for I = 2. Therefore \GiC(l,m)\ for I — 1 is much larger than for I = 2, which 
allows us to keep only the terms corresponding to I = 1 and I = 2 in the decomposition of <p 
on spherical harmonics. 
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r]{r) = lO" 3 exp(25 m 2 ) 

rip 


V = 1.26 x 10 21 Js- 1 


r)(r) = 10°exp(25 (^) 2 ) 


V = 2.7 x 10 21 Js" 1 


r](r) = 10 3 exp(25 (^) 2 ) 


P = 2.18 x 10 21 Js- 1 


r)(r) = 10 5 exp(25 (^) 2 ) 


V = 1.71 x 10 21 Js' 1 


r)(r) = 10 7 exp(25 (^) 2 ) 


V = 1.12 x 10 21 Js' 1 


V {r) = 10 9 exp(25 (^) 2 ) 


V = 2.33 x 10 20 Js' 1 


V (r) = 10 10 exp(25 (^) 2 ) 


V = 2.5 x 10 19 Js' 1 


V (r) = 10 12 exp(25 (^) 2 ) 


V = 2.5 x 10 17 Js" 1 



Table 1: Table giving V as a function of 77, where the value of 77 is artificially modified from 
the value computed in section §4 



These results indicate that the energy dissipation rate is insensitive to a change in 
the amplitude of the conductivity by several orders of magnitude (this conclusion is in 
agreement with a conjecture that Campbell made (C83)). A high conductivity increases the 
energy dissipation in a given volume, but it also tends to prevent the magnetic field from 
penetrating inside the planet. On the other hand, a lower conductivity corresponds to less 
dissipation per unit of volume, but it also allows the field to penetrate deeper inside the 
planet (and therefore increasing the volume where energy can be dissipated). 



The skin depth (for reasonable values of r\ (r)) is of order 5 = */^\ For 77 (r 



10 3 exp 



25 



(i? P ) 



and uj = 10 5 s 1 , we have S(r pn ) pa 4 x 10 7 m (we define r pn to be 

the radius of penetration, or the radius to which the magnetic field can diffuse inside the 
planet). This estimate is consistent with the numerical values of Gi(r) inside the planet (see 
figure 2) in which we find that G±(r) for r < r pn ~ 6.5 x 10 7 m is negligibly small compared 
to its value elsewhere. 

These considerations suggest that the total rate of energy dissipation is well determined 
though the location where it occurs is less well established due to the uncertainties in 77. 
Moreover, with our definition uj = u± — Q p , uj is positive outside corrotation and negative 
inside corrotation. However the ohmic dissipation rate inside the planet V only depends on 
the absolute value of uj. 



-17- 



4.3. Energy source and direct influence on the planet's orbit. 

The induced current J deduced in the previous section is due to the diffusion of a time 
dependent magnetic field. This time dependence comes from the relative motion of the 
planet's orbit and the stellar magnetosphere. Thus, the ohmic dissipation must be supplied 
by the orbital kinetic energy of the planet and the rotational energy of the star. Our stated 
goal in the introduction is to consider whether the migration of some planets may be halted 
by their magnetic coupling with their rapidly spinning magnetized host stars. In the case 
where a;* > fl p , the rotational energy of the star is transferred to the total orbital energy of 
the planet and provides a supply for the ohmic dissipation. The torque T associated with 
the ohmic dissipation is linked with the ohmic dissipation rate V and the relative angular 
velocity u according to the following equation (cf. C83, eq. 55), 

V = -lo\T\. (16) 

Since the transfer of angular momentum requires the torque associated with the ohmic 
dissipation, a similar fraction of energy is being transfered to the planet's orbit and supplied 
to the ohmic dissipation. For this purpose, we qualitatively compare the magnitude of V(t) 
with the rate of energy change needed stall the migration of a protoplanet. A detailed 
computation on the orbital evolution of the planet will be presented in Paper III. 

For illustration purpose, we first consider the power associated with the migration (P m i g ) 
of a planet with a 0.63 Jupiter mass and a 1.4 Jupiter radius toward a sun-like star. At any 
semi major axis a, the total energy of the Keplerian orbit is \E\ = GM p M s /2a. If its orbit 
decays on a characteristic planet-disk interaction time scale (r m i g ) of about 3 million years, 
the torque needed to halt the planet's migration would correspond to a power P mig such that 

P mig = \E\ ~ ~ 7.4 x lCPJs- 1 . 

Since this power is more than an order of magnitude larger than the time average value of 
V (see Table 1), it seems, therefore, not possible for the magnetic coupling to directly stall 
the planet's migration at a 0.04 AU Keplerian orbit within a few millions years, even in the 
limit of a positive uj. 

However in the model we have considered here, V oc B 2 oc m 2 a~ 6 . The power needed to 
drive the planet to migration with a specified speed is proportional to a~ 2 (these scalings are 
confirmed by numerical calculations that neglect any changes in the relative frequencies u and 
the planetary internal structure). It means that there is a semi-major axis a stop (~ 0.01AU) 
at which V and P mig are comparable. This distance is comparable to the radius of a typical 
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T Tauri stars. Note that the requirement for > f2 p also implies that the planet must 
be outside the corotation radius. This condition is satisfied only in a disk with a low gas 
accretion rate around a rapidly spinning and weakly magnetized stars. In paper II, we 
will consider such a model for the newly discovered planet around TW Hyd (Setiawan et 
al. 2008). Under these circumstances, the planet-star magnetic interaction may also be 
overwhelmed by their tidal interaction. 

5. Planetary inflation and mass loss 

In this section, we propose that ohmic dissipation in the planet's interior can indirectly 
halt its migration. The main physical mechanisms involve the heating of the planet's interior, 
its inflation and mass loss through Roche lobe overflow, and angular momentum transfer from 
the transferred material to the orbit of the planet. 

Up to now, we have computed the planet's conductivity for one particular set of param- 
eters (M p , a, etc.), and the corresponding ohmic energy dissipation inside the planet due to 
the star's magnetic field. Although, this dissipation rate for most close-in planets is generally 
too small to directly provide the power needed to halt their migration over the time scale of 
a few Myr, it can modify their internal structure. 

The ohmic dissipation is likely to increase the temperature, the ionization fraction, and 
the conductivity around the region where most of the dissipation occurs. In principle, the 
extra energy source would reduce the skin depth. However, the envelope of the young planet 
is likely to be fully convective, similar to the low-mass main sequence secondary in interacting 
binaries. Campbell (C83) suggested that the dominant diffusivity may be due to turbulence 
(Cowling 1981). In §4.21 we have already indicated that even though the skin depth may be 
affected by the magnitude of the diffusivity, the total energy dissipation rate in the planet's 
interior is not sensitively determined by the profile of r]. 

Nevertheless, the heat released by the dissipation is comparable to that associated with 
the Kelvin- Helmholtz contraction during the early stage of the planet's evolution (Boden- 
heimer, Lin, Mardling 2001 (BLM)). In the proximity of its host star, this extra energy 
source may cause a planet to inflat beyond its Hill's radius and lose mass (Gu et al. 2004). 

In the following sections (§§5-7), we adopt an idealized and self-consistent model of 
the planet's internal structure. This approach allows us to compute the conductivity of the 
planet for different sets of parameters. Considering the low dependence of the total ohmic 
dissipation V on rj, an idealized but versatile prescription is adequate for the computation of 
V and the mass loss rate (M) for different values of the important parameters (the planet's 
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mass and radius, the star's mass, luminosity, and dipolar magnetic field strength, the tilt 
between the magnetic dipole and the stellar spinning axis, and the relative orbital period). 

In §5, we show how the mass loss rate M is related to the ohmic dissipation V. In §6, we 
describe the model we used for the planet's interior and in §7 we calculate V and M for 
different sets of parameters. 

5.1. A qualitative description 

The planet receives energy, at its surface, from the star's radiation and, in the interior, 
from the ohmic dissipation. The surface heating diffuses inwards until an isothermal structure 
is established in the planet's outer envelope. But, well below the surface region, the heat 
flux is generated by the planet's Kelvin-Helmholtz contraction and ohmic dissipation and 
transported by convection. In the limit that convection is efficient, the envelope attains a 
constant entropy profile. For computational simplicity, we adopt an isothermal model near 
the surface of the planet and a polytrope model for its deep interior. 

There are two regions of interest. Very close to the host star, the ohmic dissipation 
rate is larger than that (L p = L*(R p /2a) 2 ) due to the stellar irradiation (L*) received by the 
planet. In this limit, the planet would rapidly expand beyond its Roche lobe and become 
tidally disrupted. In accordance with the results of the previous section, (in which the effect 
of V on the internal structure of the planet has been neglected), V oc cr 6 and L p oc cr 2 . 
Thus, the stellar heating dominates at larger semi major axis. In this section, we consider 
the effect of planet's inflation due to the ohmic dissipation and show that V also increases 
with the planetary radius R p at nearly the same rate as L p ((x R 2 p ). Thus, during the thermal 
expansion of the planet, the ratio of L p /V does not change. In the region where L p > V, 
the effective temperature at the planet's surface, with or without the contribution from 
the ohmic dissipation remains to be the equilibrium value T p . But, the planet's radius for 
thermal equilibrium increases with V which adds to the energy generation in the planet's 
interior (BLM). If the new equilibrium R p is larger than the planet's Roche radius, Rh, mass 
would be lost gradually through Roche overflow. 

5.2. Mass loss rate 

We now derive the equations that will allow us to calculate the mass loss rate M and 
angular momentum transfer rate as functions of V. We are in the second region where the 
ohmic dissipation is less than the radiation flux from the star (V ^ L p ), and we set the 
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Bond albedo to zero. We therefore assume that the equilibrium temperature at the surface 
of the planet is fixed by the radiation from the star T p = w ^ st a ar a2 (L star the total luminosity 
of the star, a r = 5.67 x 10~ 8 Js~ l m~ 2 T~ A ), and that the ohmic dissipation provides the 
additional energy to innate the planet. 

An irradiated short-period planet establishes an isothermal surface layer. The hot inte- 
rior continues to transport heat to this region and then radiates to infinity with a luminosity 
Lj despite the surface heating. Note that 

Li « L p = 4ira r TpRp (17) 

so that the modification to T p is negligible. The magnitude of Li is a function of R p , M p , 
T p , and the existence of the core. We have previous computed an equilibrium model for 
the parameters for several short-period planets (BLM). In the range [lO -8 Z/ lO -5 I/0], the 
numerical results of BLM can be approximated by 

log^ = A(M p ) + B{M p )\og^- + C(M p ) (log^-) ' . (18) 

Kq Lq \ Li® J 

For HD209458b (the 0.63M/ model we presented in the previous section), (A, B, C)= (3.11, 
1.01, 0.0642). BLM also determined the value of these coefficients for more massive planets 
around a solar type star (they are modified by the stellar irradiation so that they are also 
function of M*). The planet's radius R p would contract unless there is an adequate energy 
source to replenish its loss of internal energy. If the ohmic dissipation can provide such a 
source, R p = R e and Li = V in a thermal equilibrium. 

At a=0.04 AU, the Roche radius of the planet is 

R H = a ( —f . 19 

From equation (JIB]), we find that the equilibrium R e ~ Rjj if L^ ~ 1O~ 5 L which is approxi- 
mately the value of P(~ 10 21 J s _1 ) we have determined for HD209458b. During the planet 
and star's infancy, this planet would inflate to fill its Roche lobe when it has migrated to 
this location. 

Outside ~ 0.04 AU, V decreases rapidly with a. Consequently, R p reduces to the value 
which is essentially not modified by the ohmic heating. For the calculation presented in the 
previous section, we neglected the inflation of the planet. In the next section, we construct 
a self consistent model taking into account of the modification of the dissipation rate due to 
the internal structural changes. For a = 0.04 AU, the intense ohmic dissipation rate (with 
Li = V ~ 2 x 10 21 J s _1 ) modifies the planet's internal structure and inflates its radius to 
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R e ~ O.5i2 . The inflation is more severe at a < 0.04 AU because V is a rapidly decreasing 
function of a. If i? e > Rh at this location, the planet would overflow its Roche lobe and loss 
mass. For the rest of this paper, we assume that the we are in the case where the planet fills 
its Roche lobe, i.e. R p = Rh 

Two remarks are appropriate here. Firstly, in order for the Roche lobe overflow to pro- 
vide angular momentum, the actual shape of the Roche lobe should be taken into account. 
However, for computational simplicity, we adopt in this paper a spherically symmetric ap- 
proximation (refer to Gu, Bodenheimer, Lin 2003 for a detailed study (GBL)). Secondly, 
we have only considered the contribution of ohmic dissipation V to the planetary inflation. 
The tidal (gravitational) interaction between the star and the planet can also significantly 
enhance the planet's inflation in some cases. More precisely, this tidal interaction can be 
strong for small semi-major axis (the tidal effect varies as a _13//2 ) and for large radius (thus, 
the more inflated the planet, the stronger this effect becomes). 



5.3. The governing equations 

Mass loss process is initiated when R e > Rh- In this limit, a contiuous flow would be 
established in which the inflation of the envelope's drives a steady supply of gas to the Roche 
lobe region. Well inside the Roche lobe, the gravitational potential is primarily determined 

by the mass of the planet M p : <\> g = — E , but near Rh, we need to take into account of 

both the planet and the star. In a frame which corotates with the planet, the gravitational 
potential U(r) 

where r is the distance to the center of the planet. The value of Rh is determined from 
dU/dr(R H ) = 0. 

In principle, this potential introduces a complex multi-dimensional flow pattern, espe- 
cially near the Roche lobe. But the expansion of the envelope originates deep in the envelope 
where the ohmic dissipation occurs. In this region, spherical symmetry is adequate. Near 
the Roche lobe, we adopt the results obtained by GBL. For computational convenience, we 
neglect the planet's spin. 

We consider a low-velocity quasi-hydrostatic expansion of the envelope. Under this 
gravitational potential in eq. (120]) . the radial component of the hydrodynamics momentum 
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equation for a volume of gas is reduced to: 

dv 1 dP. 1 

where v is the radial velocity. The radial component of the equation of mass conservation: 
V(p(r)-u(r)) = gives r 2 p(r)v(r) = Constant, and then the mass loss rate M is constant: 

M= 4vr r 2 p{r)v{r) = Constant (22) 

or equivalently, = Then using ^£ = = &c 2 s (c 2 (r) is the sound speed), 

the momentum equation becomes: 

v 2 ) dr\2 ) ' pdr I r f 1 ' 

At a (sonic) radius r 2 near the inner Langragian point, the flow velocity becomes comparable 
to the sound speed (GBL) i.e. 

v(r 2 ) = c s (r 2 ). (24) 
where the magnitude of r 2 is the largest solution of (r 2 — rR H + j^jjf) = 0, with 

M n \ ( 2a 3 \ M„ / ' / a x 3 



^l 1 "^ l(^w + 1 J + ^i 2 U) +1 ^ 9 ' (25) 

The expansion rate is determined by the rate of ohmic energy dissipation within the 
planet. In a steady state, the energy equation reduces to 



1 d 

r 2 dr 



v 2 

r 2 p(r)v(r) ( — + h(r) + (j) g (r) 



V vol (26) 



where V vo i is the volumic ohmic energy dissipation, and <p g the gravitational potential (of 
the planet only or of both the planet and the star, depending on the location). In this 
approximation, we assume that the distribution of enthalpy h is determined by both efficient 
convective transport (in term of an abiabat) and radiative diffusion inside the planet. 

Using equation ( |22l) . we replace r 2 p(r)v(r) by ^ in eq. [23 We then can then integrate 
equation (126]) between the radius r pn (the radius at which the field can no longer penetrate 
into the planet) and r 2 so that 

M (y( r 2) - y (?>) + h(r 2 ) - h(r pn ) + 9 (r 2 ) - <j) g {r pn ) \ = J 4irr 2 V vo i(r)dr = V. (27) 
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Within an order of magnitude, A(/i) ~ — | A (0 9 ) and A(|t> 2 ) « ^ A (</> g ) (this comes from 
calculating the order of magnitude of these 3 terms using an order of magnitude for the 
temperature, for the sound speed, and for r pn ). In addition, the integrated energy equation 
is the result of an approximation as the total ohmic dissipation rate V should be the integral 
of V vo i between r pn , and R p (= Rh because we assumed that the planet fills its Roche 
lobe). However, this approximation is reasonable since h(Rjj) = h(r 2 ) (the surface region is 
approximately isothermal), w 2 (r 2 ) ~ v 2 (Rh), and </> g (r 2 ) ~ (P 9 (Rh because r 2 is very close 
to R H ). 

We now can calculate p(r 2 ) and -P(r 2 ). Equation (1221) for r = r 2 with v(r 2 ) = c s (r 2 ) = 
y^4 meters gives: 

P(r 2 ) = a p(r 2 )T(r 2 ) ( 28 ) 

where J\f a is the Avogadro constant, ks the Boltzmann constant, M.h is the hydrogen molar 
mass, /i a coefficient which depend on the ionisation rate (/x = 1 for hydrogen atoms, and 
/i = 0.5 for fully ionized hydrogen gas, and we usually choose fi close to unity). 



6. Isothermal and polytropic model 

In §4, the permeation and dissipation of the time-dependent external field is analyzed 
by neglecting any resulting changes in the planet's interior. In §5, we show that the resulting 
ohmic dissipation can substantially modify the temperature and density distribution within 
the planet. Increases in the ionization rate modify the skin depth and relocate the region 
of maximum ohmic dissipation. However, the expansion of the planet's envelope does not 
affect the rate of ohmic dissipation. In this section, we present a set of approximately self- 
consistent calculations to analyze the feedback effect of ohmic dissipation on the planet's 
internal structure. 



6.1. A Roche-lobe filling structural Model 

In principle, the structure of the planet should be solved numerically with the standard 
planetary structure equations (BLM). However, a semi analytic model based on simplifying 
assumption may provide insight on the inter-dependent relation between various physical 
parameters. Based on the BLM's numerical models, we approximate the internal structure 
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of the planet with an idealized model in which the outer region is isothermal (due to the 
stellar irradiation) and the inner region is polytropic (due to an efficient mix of entropy by 
thermal convection). In the computation of rj, we only take into account the ionization of 
the hydrogen because the internal temperature distribution is mostly determined by heat 
transfer rather than heat dissipation and the rate of V is a relatively insensitive function of 
f]. The advantage of this approximation is that its application for the self consistent analysis 
is relatively straightforward. 



The isothermal region: it extends from the surface to a transition radius r= which is to 
be determined self consistently in §6.21 In this region, the equation of state and the 
equation describing the hydrostatic equilibrium are: 

T(r) = Constant, (29) 

P(r) = a p(r) T(r), (30) 

dP_,, = GM mt (r) 

dr r 2 



(r) = -^f^p(r), (31) 



where a = and M int (r) is the planet's mass inside a sphere of radius r centered 
on the planet's center. For all practical purpose, p is sufficiently low in the isothermal 
region that we can approximate M int (r) ~ M p i anet (one can verify, a posteriori, that 
the neglected mass is less than a few percent of the total mass). To calculate P(r), we 
integrate (IBTj) using (T2"9l and (13"Tj) . We then can calculate p(r) using (|3"T|) : 

P(r) = C exp i) 
p{r) = JL C exp JL,) 

where C is an integration constant, which value is obtained by injecting T2 in the 
previous equations. 

Polytrope region: it extends from the center of the planet to r = . In this region, we use 
the following equations: 

P(r) = Kpi(r), (33) 

%r) = GM 7 (r) , (34) 
dr r 2 

dP GM mt (r) 

= ~2 ( 35 ) 

A<p g (r) = 4vr G p(r) (36) 

where <p g is the gravitational potential, and A is the laplacian (in the Poisson equation). 
Using equations ( |33l) and (1341 . equation ([35]) becomes: i^7p 7_1 (r)^(r) = — p(r)^|f 
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And after integration: <Pg{ r ) — Constant — 
We then replace <p g in the poisson equation (l36l : 



Ap 7 " 



(r) 



7 



1 



47rGp(r) 



(37) 



X 7 



For the condition appropriate in the interior of planets, the equation of state is rea- 
sonably approximated by a 7 = 2 polytrope (de Pater & Lissauer 2001). In spherical 
coordinates, the previous equation becomes: 



This equation has an analytical solution (Ogilvie & Lin 2004), and we can calculate 
p(r), P(r) and T(r) in the region described by the polytropic equation of state: 



In principle, the transition between the two regions is determined by the onset of con- 
vection. In the construction of hydrostatic equilibrium structure models (to be presented 
in Paper II), we will indeed use that condition to determine its photospheric radius. Qual- 
itatively, we expect the transition radius which separates the two regions, r= to be larger 
than r pn , because only in the region interior to r= do we expect the temperature, ionization 
fraction, and conductivity to be sufficiently large to halt the penetration of the field. In a 
hydrostatic equilibrium, the actual value of r= is determined by the ratio of the ohmic dis- 
sipation rate in the convective region to the sum of the ohmic dissipation rate in the entire 
planet's interior and the stellar irradiative flux on the planet's surface. A set of fully self- 
consistent solution requires the matching of the ohmic dissipation rate to be expected from 
the planetary structure and that which determines its density and temperature distribution 
(see paper II). 

In the present context, we are considering the situation in which the planet's radius 
is constrained by its Roche lobe and the density and temperature of the outer boundary 
is determined by equation fT28l . In this configuration, heat is also transported by advection 
which modifies the location of r = . Moreover, the density ratio between the planet's center 
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6.2. 



Transition between the two models 
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and the outer boundary is much larger than the temperature ratio. Therefore, the polytropic 
region cannot fill the entire interior region. Since the pressure scale height on the planet's 
surface is much smaller than its radius, the isothermal region also cannot occur in the entire 
planet's interior while containing all of its mass. Instead, the planet's interior adjusts to 
attain a balance between the requirement of mass loading and constraints set by hydrostatic 
equilibrium for appropriate equations of state. 

In order to construct such an equilibrium model, we now determine po an d k at r = where 
the transition between the two regions occur. There are three equations that constrain 

these parameters. Tisothermal y=) = T p olytro P ic{j'=)i Pisothermalij '=) = Ppolytropicy '=); and the 

total mass is constant. The first two conditions also imply Pi SO thermai(r = ) = p po iytro P ic(r=). 
Therefore, we solve the following equations for p , k, and r=: 



By assuming an isothermal structure in the outer envelope, we have neglected an out- 
ward heat flux. This approximation is only adequate if the dissipation rate is above that 
which is need to inflat R p to the planet's Roche radius. If this condition is not satisfied, the 
planet's radius would attain equilibrium values for which the surface cooling is balanced by 
the Ohmic dissipation and stellar irradiation. We will construct, in Paper II, the equivalent 
of equation 21 (for a 0.63 M & planet) which takes into account the effect of ohmic dissipation 
in the planetary interior. 

Whereas the temperature on the planet's surface is determined by the stellar irradiation, 
the density at r 2 = Rh is determined by the magnitude of M (through equation which 
in term is determined by the rate of ohmic energy dissipation V (see §7]). 

For very large values of V, a set of fully self-consistent solutions also modifies the 
temperature at the disk surface as well as the thermal content of the outflowing gas. However, 
provided V is small compared with the stellar irradiative flux, a transition for convective 
stability occurs near r = . 




r o a 47rr 2 p P oiytropic(r)(ir + J a 4:Trr 2 p isothcvmal (r)dr = M p . 



(40) 



6.3. Calculation of the magnetic diffusivity 



With these internal structure specified, we consider Saha's equation for the hydrogen 
atoms which gives the ionization fraction x (Kippenhahn & Weigertal, pages 107-111): 



where the ionization energy of hydrogen is E = 13.QeV. We also neglect here the radiation 
pressure as we write P gas ( r ) = P( r )- 

If the ionization fraction x is small, x 2 ~ Kh (this is typically the case in the region 
where the ohmic dissipation occurs). We use a = with v e = N n 10~ 19 )^ . 

The electric conductivity we would obtain does not take into account higher ionization 
states or the ionization of elements other than hydrogen atoms. We then use for the following 
calculations an electric conductivity that is 10 times higher than that we would obtain with 
the Saha equation (eq. f|4T]) ) for the hydrogen atom only. We saw in table 1 that the ohmic 
dissipation rate was quite insensitive to the magnetic diffusivity i](r) = (fi cr(r)), and we 
verified that this is also the case with the internal model we used for the planet in sections 
5 to 8 (for example, in this model, a uniform change in the magnetic diffusivity by a factor 

10 changes V and M by less than 20%, and a uniform change in the magnetic diffusivity by 

a factor 100 changes V and M by less than 40%). 

We then obtain the following expression for the magnetic diffusivity inside the planet: 



where T(r) and P(r) are the temperature and pressure of the isothermal or polytropic region, 
depending on the radius r. 

7. Ohmic dissipation rate and the mass loss rate for different sets of 



With the above idealized prescription for the planet's internal structure, we now calcul- 
tate self-consistently the ohmic dissipation rate V inside the planet as well as the mass loss 




(41) 




(42) 



parameters 



rate M. 
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7.1. Parameters involved in the calculation 

The model parameters involved in the calculation of the ohmic dissipation rate inside 
the planet are: 1) the planet's mass M p , 2) semi-major axis a, 3) the relative angular velocity 
uj (the angular velocity of the field seen in a frame centered on the star and rotating with 
the planet), 4) the strength of the star's magnetic dipole moment m, and 5) the angle a 
between the spin axis of the star and the star's magnetic dipole. 

We use the isothermal and polytropic prescription described in the previous section to 
model the planet's internal structure and calculate the conductivity profile inside the planet. 
To do so, we also need to specify 6) the mass of the star M*, and 7) the star's luminosity L*. 

7.2. Methodology 

The construction of a self-consistent model requires a loop of retroaction involving the 
determination of the internal structure of the planet and that of the ohmic energy dissipa- 
tion. For a specified internal structure of the planet, one can compute (following §3) the 
conductivity profile and then the total ohmic dissipation rate V . However, this energy dissi- 
pated inside the planet corresponds to an input of heat, which triggers an adjustment in the 
planet's internal parameters. Due to the efficient convection inside the planet, we assume 
that the adjustment of the internal parameters to this external heating is quick. We consider 
that the characteristic time scale for the planet to evolve from one equilibrium state to an- 
other is small compared to the variation time scale of the seven parameters mentioned in the 
previous paragraph. Therefore, we do not need to follow the planet's dynamical evolution 
at all times. Instead, we can take a series of "snap shots" of the planet in its equilibrium 
state for different set of parameters. 

Because of this feedback loop between the ohmic dissipation rate and the planet's inter- 
nal parameters, we use an iterative method. For any chosen set of parameters, we start from 
an estimate for the ohmic dissipation rate Vo and internal structure T (r), Po(r), and po(r) 
corresponding to a magnetic diffusivity profile r?o( r ) (I n our parametric analyses, we typi- 
cally make small incremental changes in the model parameters from those for which we have 
already obtained equilibrium values). We then compute the new internal structure Ti(r), 
-Pi(r), and pi(r) associated with Vq. This enables us to compute the corresponding magnetic 
diffusivity rji(r). Finally, we use rji (r) to calculate the corresponding ohmic dissipation rate 

V\ and mass loss rate M\. This process is iterated until convergence of V, M, and of the 
internal parameters. Moreover, for some specific set of parameters M p , a, u, m, sin(a), M*, 
L*), we have started the iterative process from two different initial states in order to verify 
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that they both converge to the same solution. Therefore, the iterative process does converge 
to a unique solution. 

We consider the following fiducial model in which the mass of the planet and semi-major 
axis corresponds to HD 209458 b, and in which the other parameters are reasonable ones for 
the type of systems considered. An estimate of the order of magnitude for the strength of 
the magnetic dipole can be found in Johns-Krull 2007. 

Mass of the planet: M p = 0.63 Mj = 1.26 x 10 27 kg, 
Semi-major axis: a = 0.04 AU = 6 x 10 9 m, 
Relative angular velocity: uj = 10~ 5 s _1 , 
Star's magnetic dipole: m = 4 x 10 34 A m 2 , 
Value of the tilt of the magnetic dipole: sin (a) = 1, 
Mass of the star: M, = M = 2 x 10 30 kg, 
Luminosity of the star: L t = 1.5 L = 5.7 x 10 26 W. 



7.3. Computation of V and M, plots, and mathematical relations 

We present seven groups of plots (figures 3, 4, 5, 6, 7, 8, 9), one group for each param- 
eter mentioned just above. For each group, we vary one parameter (x-axis), while keeping 
the others at the reference values mentioned above. On the y-axis, we plotted the ohmic 

dissipation rate V, mass loss rate M, and characteristic time scale tm = —■ Note that the 

M 

magnitude of tm ~ for a Jupiter mass planet is about lMyr. In addition, the mass loss rate 
determined here is many orders of magnitude larger than that due to photo evaporation. 
Only with such large mass loss rate, can we compensate for the angular momentum transfer 
due to the planet-disk and planet-star tidal interaction. 

We emphasize once again that in the construction of these models, we assume that there 
is adequate energy dissipation to inflate the planet with R e > Rh- In later papers that use 
the roche-filling model, we verify that V(Rh) > L^Rh) before these results are applied. If 
this condition is not satisfied, the planet would not fill its Roche lobe and not lose mass. 

From each group of plots, we obtain V and M as a function of the parameter that is 
being varied (all the others are kept constant at the value of the fiducial model given above). 
These functions are given in the table below (table [2]). 
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Ohmic dissipation rate V and mass loss rate M 

TW 



Varying Parameter 



V x = 3.3 x 10 



21 



Mi= 1.2 x 10 



13 



( p | w 

V 0.63 Mj J vv 
\ 2 4 



0.25 Mj < M p < 1.7 Mj 



0.63 Mj 



kg s 1 



P 2 = 3.3xKF {^wY W 
M 2 =1.2xl0 13 (^^kgs' 1 



0.015 AU < a < 0.08 AU 



V 3 = 3.5 x 10 21 j - 1.5 x 10 20 W 



7 x 10~ 6 s" 1 < U < 7.3 x 10~ 5 s- 1 



M 3 = 1.4 x 10 



13 



10" 



1.3 x 10 12 kg s 



-i 



W 



V 4 = 3.3 x 10 21 (^) 
M 4 =1.2x 10 13 (^^kgs' 1 



6 x 10 33 Am 2 < m < 4 x 10 34 Am 2 



V 5 = 3.3 x 10 21 sin 2 - 17 (a) W 
M 5 = 1.2 x 10 13 sin 2 - 28 (a)kg s' 1 



0.3 < sin(a) < 1 



n = 3.3 x 10 21 (^j \Y 
M 6 =1.2x 10 13 ^j'^kgs' 1 



0.5 M < M* < 1.5 M P 



•^v 3.3 x 10 21 (nrfeP^ 
M 7 =1.2x 10 13 °' 8 kgs 

5.9 



-1 



V 7 = 7.5 x 10 



19 



1.5 L 



M 7 = 1.25 x 10 



ii 



1.5 L P 



w 

kg s _1 



0.5 L Q < L*< 2.6 L Q 



2.6 L Q <L+< 5 L 



Table 2: Table giving P and M as a function of the parameter that is being varied. 
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Fig. 3. — Ohmic dissipation rate, mass loss rate, and time scale for different planetary masses 

7.4. Model parameter dependence 

1. Mass of the planet M p . The total ohmic dissipation is a volumic integral over the 
entire region where dissipation occurs. Therefore, one might expect V to be propor- 
tional to the volume. Since the planet fills its Roche lobe, the volume is detemined by 
the mass (cf. eq. fT9l) . However, M p also gives a constraint on the volumic mass at 
the center of the planet (cf. equation SOD , which makes V mostly proportional to M 2 . 
There is also a minor correction due to Vs weak dependence on rj which depends on 
M p through the calculation of the internal parameters T, P, p (eqs. [321 and [391) . 

2. Semi major axis a. In the model we adopted in §4, or more generally, in a model that 
would not take into account the planet's internal adjustment to the ohmic dissipation 
(especially in a model in which the radius of the planet is independent of the semi-major 
axis), the ohmic dissipation rate would be related to the semi-major axis according to 
the following law: V oc B 2 oc a -6 . In the self-consistent model we adopted here and 
in the limit where the planet fills its Roche lobe (R p = Rh), the radius of the planet 
varies with the semi-major axis. For example, when a planet moves closer to its host 
star, its Roche radius decreases linearly with the radius (cf. eq. [19]) and, therefore, 
V increases less quickly than if the planet kept the same radius. Again, V also has a 
weak dependence on 77 which is depends on the semi-major axis through the planet's 
surface temperature and through the dependence of on a. From these arguments, 
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Fig. 4. — Ohmic dissipation rate, mass loss rate, and time scale for different semi-major 
axies 

we thus expect V the exponent in Vi to be greater than -6 and less than -3. 

3. Relative angular velocity uj. In equation (??, the multiplicative constant in front of 
the volumic integral comes comes from the induction by the time dependent stellar field 
(cf. eq. (JH]) and eq. (fTTj) ) and gives to V a dependence on uj 2 . However, in addition, 
uj intervenes inside the volumic integral in eq. fT2l through the dependence of Gi(r) on 
the 5" 1 , with 5 = a/J, skin depth. Therefore, the volumic integral is proportional to 
uj^ 1 and V is proportional to uj. 

4. Magnetic dipole m and tilt of the magnetic dipole a. Without any adjust- 
ment of the planet's interior to the ohmic dissipation, we would expect P4 and V$ to 
vary respectively in m 2 and sin 2 (a). The fact that both exponents that have been 
computed numerically are slightly larger than 2 means, in the self-consistent model 
we adopted here, that the adjustment of the planet's interior tends to have a small 
positive retroaction on the amount of enegy that is deposited inside the planet through 
ohmic dissipation. 

5. Mass of the star M*. The mass of the star intervenes in the computation of the 

— 1/3 

Roche radius (Rh oc M* ) . V being a volumic integral, one would, therefore, expect 
it to vary proportionally to M~ l . However, M* also intervenes in the computation of 
r 2 (r 2 is the sonic point, or the largest solution of (r 2 — rR H + jTfrr-) = 0, see. eq 
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Fig. 5. — Ohmic dissipation rate, mass loss rate, and time scale for different relative angular 
velocities 

1251) . which brings a correction to the dependence of V on M s . Indeed, r 2 is used to 
compute P(r 2 ) and p(r 2 ) (cf. eq. [28]) which are the boundary conditions we adopted 
to calculate the pressure and volumic mass in the isothermal region (cf. eq. [321) . As a 
side note, this dependence of V on r 2 could also affect the dependence of V on a and 

6. Stellar total Luminosity L*. From figure 9, one can see that L* = 10 27 W cor- 
reponds to a minimum for V and that V varies slowly for < 10 27 W and much 
faster for > I0 27 W . In the model we adopted here, the stellar total luminosity 
fixes the planet's equilibrium surface temperature, which is also the temperature of 
the isothermal region (for = 10 27 W, T p ~ 1767K.). It in turns determines the 
temperature profile inside the planet (the surface temperature is used as a boundary 
condition) and influences the internal structure and magnetic diffusivity profile inside 
the planet rj(r). T p varies proportionally to L l J A (for constant semi-major axis) and, 
therefore, r\ is roughly proportional to exp(^x). 

We ploted the integrand of the ohmic dissipation (< V vo i um i C > r 2 in eq. [13]) as well as 
the magnetic diffusivity (see figures 10 and 11) for L t = 2 x 10 26 W, 7 x 10 26 W, 10 27 
W, 1.2 x 10 27 W, 1.5 x 10 27 W, and 2 x 10 27 W respectively. One can notice two parts 
corresponding to the isothermal and the polytropic region. In the isothermal region, 
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Fig. 6. — Ohmic dissipation rate, mass loss rate, and time scale for different stellar magnetic 
moment 

the energy integrand decreases slowly from the surface to the center of the planet. 
The transition between the two regions corresponds to a sharp increase in conductivity 
and, therefore, a quick increase in the ohmic dissipation. This accounts for the sharp 
increase in the integrand (around 1.5 x 10 8 meters), and most of the remaining mag- 
netic energy is dissipated in this region. Moreover, an increase in the stellar luminosity 
results in a decrease in the magnetic diffusivity as well as a deeper penetration and a 
deeper transition between the isothermal and the polytropic region. 

When one increases starting from low values (L+ = 2 x 10 26 W), the energy integrand 
also increases in the isothermal region. Nevertheless, the extent and amplitude of the 
sharp increase at the transition between the isothermal and the polytropic region are 
also reduced. Therefore, these two effects compensate each other, and for low stellar 
luminosity (e.g. L* between 2 x 10 26 W and 10 27 W, the total ohmic dissipation in the 
planet increases slowly with the stellar luminosity). 

On the other hand, for higher values of the stellar luminosity, the penetration depth as 
well as the transition depth saturates (when the coupling term in eq. (jBJ) between 

the real and imaginary parts of Gi(r) becomes comparable to the other term ( "t 1 ^ )). 
One can see that an increase in the stellar luminosity results only in changes in the 
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Fig. 7. — Ohmic dissipation rate, mass loss rate, and time scale for different tilt of the stellar 
magnetic dipole with regard to the stellar spin axis 



energy integrand that would increase the total ohmic dissipation. This result in a much 
sharper increase of the ohmic dissipation with the stellar luminosity. 

7. Attempts of generalized function. We consider the generalized expression of the 
ohmic dissipation rate and mass loss rate, in the case where the variables are separable: 



V Qen = 3.3 x 10 21 W , 
9 v 0.63 M 



r_?4_y 16 (^-v \( 

\0.63MjJ \0MAUJ [\ 

2.18 



m 



\LU\ 
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sin a) 



v2.17 



lO^s- 1 
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(43) 



M 
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1.2 x 10 13 kg s" 1 



0.63 Mi 
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\U)\ 



m 



4 x 10 34 A m z 



2.3 



(sin a) 



2.28 



M ( 



lO^s- 1 

0.5 



-0.1 

-0.8 



1.5 LQ 



(44) 



The previous formulas have been written for the first interval of Vj in table 2, but 
one can write the corresponding formulas for the second interval of by using the corre- 
sponding expression of V-j. We compute, using the iterative procedure described above 
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Fig. 8. — Ohmic dissipation rate, mass loss rate, and time scale for different stellar masses 



(§7.2), the ohmic dissipation V and mass loss rate J\4 for sets of parameters in which 
more than 2 parameters are different from the fiducial parameters. We then compare 
these values with the value of V gen and J\A gen for these sets of parameters. This test en- 
ables us to determine if the hypothesis of separation of variable is accurate or not. The 
first set of parameters that we consider is M p = l.hMj and a = 0.03AU, all the other 
parameters being kept equal to the fiducial parameters. We get: V = 4.2 x 10 22 W 
and jii = 1.5 x f0 14 kg s~\ However, using the formula of the ohmic dissipation 
rate and mass loss rate with the approximation of separation of variables, we get: 

Vg en = 7 x 10 22 W and M = 2.9 x 10 14 kg s _1 . 

We now consider a set of parameters in which all parameters are taken differnet from the 
fiducial value: M p = 1.5M,-, a = 0.03AU, w = 2.9x lCT 5 ^ 1 , m = 2 x 10 34 Am 2 , M* = 
1.5M Q , L* = O.8L (this value of u corresponds, for example to a system in which the 
planet is at keplerian augular velocity and the star has a period of 4 days). We get 

V — 1.7 x 10 22 W and M = 7 x 10 13 kg s~\ and using the formula of the ohmic dissi- 
pation rate and mass loss rate with the approximation of separation of variables, we 

get V gen = 3.5 x 10 22 W and M = 1.6 x 10 14 kg s" 1 . 



From these two tests, we deduce that the approximation of separation of variables gives 
a reasonable order of magnitude, but nevertheless, does not seem to be accurate. This 
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Fig. 9. — Ohmic dissipation rate, mass loss rate, and time scale for different stellar total flux 

results means that the exponents in the functions given in table 2 can themselves be 
a function of the parameters, and the generalized functions that describe the ohmic 
dissipation and mass loss rate as a function of all parameters can be fairly compli- 
cated. Nevertheless, even if a generalized formula is still unknown, for a given set of 
parameters, one could still compute the mass loss rate and ohmic dissipation using the 
procedure described in §§3-7. 

7.5. Mass loss and migration stalls 



From equation (I44p . we find 




The same mass loss provides angular momentum to the planet. We neglect any varia- 
tions in excentricity and assume that all the mass is accreted into the star. Using (GBL eq. 



-38- 




Fig. 10. — Integrand of the ohmic dissipation (log scale) for different values of the stellar 
luminosity 

96), we link the mass loss rate to a rate of change of semi- major axis, 

- = " 2 # <«) 

a Mp v ' 



Thus, T a = \a\/a = — 2r, 



m ■ 



These relations indicate that, within ~ 0.04 AU, the Ohmic dissipation within the 
planet may indeed generate sufficient energy to inflat their radius beyond their Roche lobe. 
The resulting mass transfer not only reduces the planets' mass but also stall their orbital 
migration. The impact of this process on the mass-period distribution of gas giants will be 
discussed in paper III. 



8. Summary 

In this paper, we applied a model described by Campbell (in the context of interacting 
binary stars) to the situation of a planet in a protoplanetary disk interacting with the stellar 
periodic magnetic field. In §3, we showed that with a well determined electrical conductivity 
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Fig. 11. — Magnetic cliff usivity (log scale) for different values of the stellar luminosity 



profile inside the planet as well as the characteristic parameters of the system (such as the 
stellar magnetic field strength and angular velocity spin, planet radius and semi-major axis), 
one can compute the total ohmic dissipation rate V(t) inside the planet and its average value 
over one synodic period. This dissipation rate gives a good estimate of the strength of the 
Lorrentz torque exerted on the planet due to the interaction between the stellar magnetic 
field and the induced current inside the planet. When the planet is outside corrotation, 
this torque will provide angular momentum to the planet from the star and slow down the 
planet's migration. In §4, we computed V for one specific set of parameters (R p = 0.63 Rj, 
a = 0.04 AU), and also showed that the conductivity profile (all the other parameters being 
kept constant) had some influence on the location of maximum dissipation, but fairly little 
influence on the total dissipate rate V. We noted that this value of V seemed too low to 
directly provide an adequate rate of angular momentum transfert to the planet to stop its 
migration toward the host star. However, this energy input can inflat the planet's enveloppe 

and trigger mass loss M through Roche lobe overflow. The mass that overflows toward the 
central star provides angular momentum to the planet (GBL). In order to estimate this mass 
loss rate, we first linked the ohmic dissipation rate to the mass loss rate (§5). Then we used 
an isotermal-polytropic model to describe the adjustment of the planet's interior to the heat 

deposited through ohmic dissipation (§6). Finally, we computed V and M at equilibrium for 
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several set of parameters (§7). A detailed calculation on the orbital evolution of the planet 
due to this process will be presented in paper III. 



We thank Drs. P. Garaud, J.E. Pringle and F. Rasio for constructive discussions. This 
work is supported by NASA (NAGS5-11779, NNG04G-191G, NNG06-GH45G, NNX07AL13G, 
HST-AR-11267), JPL (1270927), and NSF(AST-0507424). 



A. Perfect conductor moving relative to a magnetic field 



Let's consider the flux $ of the magnetic field B accross a surface S(t) that changes 
with time or moves in space. One can show that 



def 1 
dt dt 



B(r,t + dt)dS- / B(r,i)dS 



<9 B „ / 

— VA(uAB) 



dS 



with a = S(t + dt), and b = S(t). Therefore, using the mhd induction equation 



<9B 

~dt 



V A (v A B) - V A ( V A B), 



we get: 



d ^ 1 f T / \ 71 



(Al) 



(A2) 



(A3) 



which tends to zero when the electric conductivity is large (the previous integral is a closed 
integral along a close curve). Therefore, the magnetic field's flux will be constant if o is large 
enough that the second term in the right hand side is negligible. In such a case, the field 
lines will move with the body and appear to be 'frozen.' 



B. Set of linear equations 

We give below the linear set of equation that we solved for //J" 1 , v\, z/f 1 , a±, 
a 2 ,a 3 ,a 4 ), Pi,f3 2 ), and (//|, /i^ 2 , v\, u^ 2 , 71, 72, 73, 74). The values of Gi(r) considered 

are for r = R p the radius of the planet. 
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(fi 2 + fi 2 2 )Re{G 2 (R P )) + - v 2 2 )Im{G2 (Rp)) + 2 72 ^12 ^/ff = 
(-v 2 + v 2 2 )Re{G2 {Rp)) + (~l4 + IH 2 )Im{G2 {Rp)) + 2^12^ = ~^R 2 p 12^ 
{-v 2 - v 2 2 )Re{G2 {R P )) + {-A - ^ 2 )Im{G2 {R P )) + 2^12^ = -|^^12^/ff 
^ (-^1 + th 2 )MG2 {R P )) + {A ~ V2 2 )Im{G2 {R P )) + 2 7 3^12 ^/ff = 



' fi° 2 Re{G 2 {R p )) - v° 2 Im{G 2 {R p )) - fo^r^Jf = 
-v° 2 Re(G 2 (R p )) - 3lm{G 2 {R p )) - A^^/f = yf^*^^ 
fi° 2 Re{G 2 {R P )) - v\lm{G 2 (Rp)) +lh^y/f = 
-v»Re(G 2 {Rp)) - /i°/m(G 2 + = yff^3R 2 p 



C. Coefficients intervening in the expression of the ohmic dissipation rate 
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D. Conductivity and resistance of short-period extrasolar planets 

We now calculate the conductivity and resistance (the reciprocal of conductivity) of 
short-period extrasolar planets. The conductivity of the planet is determined by the density 
of charged particles. We consider separately the contributiond from the collisional ionization 
within the planet's interior and from the photoionization near its surface. 



D.l. Ionization of the planet's interior 

The cores and the inner envelopes of Jovian planets are mostly ionized, they are shielded 
by cool, mostly neutral gaseous envelopes, where the ionization is dominated by elements 
with low ionization potentials, such as sodium and potassium. 

The planetary model used in our calculation of ionization fraction is model A3 presented 
by Bodenheimer et al. (2001). It is a spherically symmetric model for the short-period 
planet around HD 209458. The following parameters are assumed: the planetary mass is 
0.63 Jupiter masses (Mj); the equlibrium temperature at the surface of the planet due 
to irradiation from the star is T s = 1360 K; there is a solid core with a density p c = 
5.5 10 3 kg m~ 3 and a mass 0.139 Mj ( = 44M e = 0.22M p ) in the center. An energy source, 
uniformly distributed through the gaseous part of the planet, with an energy input rate 
Ed = 8.5 x 10 19 J s _1 , is also imposed to take into account the effect of tidal dissipation of 
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energy. Those model parameters result in an asymptotic radius of IAIRj at t = 4.5 Gyr, 
which is consistent with the photometric occultation observations of HD 209458 (Henry et 
al. 2000b; Charbonneau et al. 2000) 

The cores and the inner envelope of Jovian planets are mostly ionized by the pressure 
ionization effect, where the Saha equation breaks down. We have to resort to various equa- 
tions of state (EOS) tables, including the equation of state tables for hydrogen and helium 
by Saumon et al. (1995). These tables are calculated using the free energy minimization 
methods, with a careful study of the nonideal interactions. They cover temperatures in the 
range 2.10 < logT(K) < 7.06 and pressure in the range 5 < logP(Nm -2 ) < 20. The 
calculations based on which these tables were constructed also include the treatments of 
partial dissociation and ionization caused by both pressure and temperature effects. Given 
the internal structure data of p, T and P for model A3, we use those EOS tables to calculate 
electron number density distribution for the inner part of the planet. In this approximation, 
we bear in mind that the ionization fraction calculated from the free-energy minimization 
are of limited accuracy. Moreover, in the interpolation regions of both the H and He equa- 
tion of state, the data have very little physical basis but are reasonably well behaved by 
construction. 

In the envelope, hydrogen and helium are mostly neutral, and free electrons are ex- 
clusively provided by thermal ionization of elements with low ionization potentials. Among 
them, potassium and sodium have highest concentrations with a relative abundance loglN^/Nx) ~ 
—6.88 and log(N^ a /NH) — —5.67 for the solar composition. The lowest ionization poten- 
tials for these two elements are xk = 4.339 eV and XNa = 5.138 eV. We identify these two 
elements to be the sources of most of the free electrons in the planetary envelope. 

We solve the Saha equations for ionization of Na & K jointly (cf. Allen 1955) 

N 1 3 

10 6 x N e = -XNa® - t: log 6 + 20.9388 (Dl) 

iV iVa 1 

and 

Xf 1 3 

10 6 x N e = -xk& - - log 6 + 20.9388 (D2) 

where O = 5040K/T, Nj^ a and N^ a are, respectively, the singly ionized and neutral number 
density of sodium (SI units), and are, respectively, the singly ionized and neutral 
number density of potassium, and N e is the electron density. 
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D.2. Conductivity and resistivity 

Using the electron number density profile, we then calculate the conductivity using 
the formulas given by Fejer (1965). Three kinds of conductivity are of interests here. The 
conductivity <7o, which determines the current parallel to the magnetic lines of force, and 
which would exist for all direction in the absence of the magnetic field, is given by 

ffQ = ^ JVef! (D3) 

B \h>i v e ) m e v e 

where e is the electron charge, to e = — f^- and Ui = ^ are the gyrofrequencies of electron 
and ion, respectively, while u e and z/« are the collisional frequencies associated with the 
momentum transfer of electrons and ions. 

The Pedersen conductivity, which determines current parallel to the electric field, is 
given by 

_ N^e ( _ u e u e \ ^ o-q mA) 

° V ~ B \uf + w ? vl + u>l) ~ 1 + (u e /u e y [UV 
The Hall conductivity, which determines the current perpendicular to both the electric and 
magnetic fields, is given by 

N eG ( UJ 2 £ ~ Voi^e/Ve) m _v 

h B \vl + ul uf + uf) l + {ujv e f 1 ] 

In the limit of low ionization fraction, v e is closely related to the mean collisional fre- 
quencies of the electrons with molecules of the neutral gas such that (cf. Draine et al. 
1983) 

=N n <av > e _„^ N n x 1(T 19 (j^j V2 • (D6) 

where N n is the number density of neutral particles, < av > e _„ is the average product of 
collisional cross section and the relative speed between electrons and neutral particles. In the 
other limit of a completely ionized gas, we take into account both electron-ion and electron- 
electron encounters. The collisional frequency of the electrons is given by (cf. Sturrock 
1994) 

v f) ~ 10 8 -°iV e T- 3 / 2 (D7) 
In the intermediate range, we use v e = maxfi/l 1 ', ui 2 ^). 

The resistance of the planet cannot be specified exactly because an unknown shape factor 
is involved. Following Dermott (1970), we use the following expression for the resistance: 

^4/"4t' < D8 > 

r 2 p Jo o-p(r) 
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where fn is a order-of-unity parameter for the geometry. In the aligned geometry we consider 
here, the planet's resistance comes from Pedersen resistivity (conductivity). The core of the 
planet is assumed to be a perfect conductor. From Eq. ID8[ we can see that it is the cold, 
mostly neutral envelope that gives rise to most of the resistance. 
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